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Abstract: We prove a large deviations principle for the largest eigenvalue 

of a class of biorthogonal and multiple orthogonal polynomial ensembles that 
includes a matrix model of Lueck, Sommers and Zirnbauer for disordered bosons 
and Angelesco ensembles. Moreover we consider matrix ensembles in mesoscopic 
physics. 
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1. Introduction 

Muttalib introduced in [18] a new model of random matrices with two-body interactions. 
This approach was motivated by physics, where the usage of a standard random matrix ansatz 
for modelling the behaviour of disordered conductors leads to a small deviation from observed 
results. Although it was known using methods of perturbation theory that a small correction 
would be necessary for exact results (see mi), it was not known how this model could be 
solved mathematically. Muttalib started at this point and contributed a solvable random 
matrix model that considered a small correction of the classical one-body model needed for 
e.g. metallic conductors. Muttalib’s model has the probability density function 

n 

q n (xi,...,x n ) = Cn \xi - XjWxf - x e j\ Y\e~ v{Xi) l S "(^i, • • • ,x n ) , (1.1) 

l<i< 7 <n 2 — 1 

where 6 is a fixed positive number, the Xj, 1 < j < n, are the eigenvalues of a n x n matrix 
X, V is a weight function specifying the concrete model. The fraction c n normalises the 
density and E is the domain of the eigenvalues. Since density (II.ip can be rewritten with the 
help of biorthogonal polynomials, Borodin [4] introduced the term biorthogonal ensembles for 
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ensembles with density fll.lj) . Since one has a repulsion between the different eigenvalues Xi 
and a repulsion between the xf’s, the point process described in (II. ip is sometimes called a 
point process with two-particle interactions. Moreover these point processes are determinantal 
point processes (see (lj Definition 4.2.11]). 

This paper is motivated by a special biorthogonal ensemble that arises from a random matrix 
model for disordered bosons that was proposed by Lueck, Sommers, and Zirnbauer in [16]. The 
model amounts to the product of a Wishart matrix and the fundamental matrix of the standard 
symplectic form, and on the level of eigenvalues, interpreted as characteristic frequencies of 
disordered quasi-particles, one obtains the joint density 

1 n 

q n , a ( x i,...,x n ) = —- JJ \xi - Xj\ \x 2 i - x]\ JJx“e _r:Ei l[ 0iO o)™(^i, • • • , x n) (1.2) 

n ’“’ r 1 <i<j<n i= 1 

with a G N U {0} and r > 0. In [16] it was shown that the correlation functions of the 
frequencies in the bulk of the spectrum are in the Gaussian Unitary Ensemble universality 
class, yet a novel scaling behaviour is found at the low frequency end of the spectrum. Other 
applications of biorthogonal ensembles to physics are discussed in [19] motivated by matrix 
models for Chern-Simons theory. For any sequence of random numbers xi,...,x n consider 
L n := the empirical distribution or empirical measure of these values (a random 

probability measure on R). Define the mean empirical measure L n = E L n by the relation 
(L n , f) = E (L n , f) for all continuous and bounded functions / : R —* R. Now in (11.211 we 
consider r = n and a 6 N U {0} be fixed. This is the case when the Gaussian random variables 
in the random matrix model in [£6] are chosen to be random variables with zero mean and 
variance 1/n . One result in [IB] is that with r = n the sequence (L n ) n converges weakly to a 
probability measure on R with Lebesgue density 

(?»(*) := 4-(i/6)-‘/3( (1 + yi - wjwy/ 3 - (i - yi-iv^w), (1.3) 

for 0 < t < b := 3\/3. Hence one would expect that the largest frequency 

* * n 
x := x n := rnaxx,- 

3 =1 

converges almost surely to the right endpoint 3\/3 of the support of g^. One aim of this paper 
is to complement the Lueck-Sommers-Zirnbauer results by a large deviations principle for the 
largest characteristic frequency x*. 

Recently, in [5] it was shown that the biorthogonal polynomials associated to models (11.11) 
satisfy a recurrence relation and a Christoffel-Darboux formula if 6 6 Q. Moreover the authors 
express the equilibrium measure associated to this model. 

Actually, we will study large deviations principles for largest eigenvalues in a broader frame¬ 
work that encompasses not only ensembles like (11.11) . but also takes care of weight functions w n , 
depending on n, and determinantal parts like n*<y \ x i ~ x j\ l3 ( x i ~ x j) f° r an y /3 > 0. Moreover 
we consider random matrix ensembles in mesoscopic physics subsuming matrix versions of all 
classical symmetric spaces. We will also obtain large deviations results for multiple orthogonal 
polynomial ensembles (see EE]). 
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In0, Ben Arous, Dembo and Guionnet proved a large deviations result for the GOE ensemble 
where the joint probability density of the eigenvalues is given by 

n 

q n (xi, ■ ■ ■ ,x n ) = c n \xi - Xj\Y\e~^ x i. (1.4) 

l<i<7<n z=l 

ffere the sequence of largest eigenvalues (x* n ) n fulhls a LDP in M with speed n and good rate 
function 

/* y/(t/ 2) 2 - 1 dt : x>2 

oo : x < 2 

Recall that a family of probability measures (/i £ ) e> o on a topological space X is said to obey 
a large deviations principle (LDP) with speed e^ 1 and good rate function / : X —> [0, oo] if 
I is lower semi-continuous and has compact level sets Nl := {x G X : I(x) < L}, for every 
L G [0, oo), and 

lim inf slog fJ, e (G) > — inf I(x) 

£—^0 X^G 

for every open G Cl and 

lim sup £ log n £ (A) < — inf I(x) 

e —^0 

for every closed A C X. The result in [2] has been generalized in [I] to joint densities of the 
form 

q n (x 1} ... ,x n ) = c n \xi - xjf e~ n ^= lV( - Xi \ (1.5) 

for (3 > 0, partition function c“ x and continuous weight functions V which satisfies 
lim infixi-^oo V(x)/(j3' log |a:|) > 1 for some j3' > 1 with /3' > /3 (see Theorem 2.6.6 in jl]). 

Recently, in [10J the author proves a large deviations principle for the largeste eigenvalue of 
Wigner matrices without Gaussian tails, namely such that the distrution tails P(|-X’ lil | > t) and 
P(|X 1)2 | > t) behave like e~ bt " and e _at “ respectively for some a, b G (0, cxd) and a G (0, 2). The 
large deviations principle is of speed n Q//2 and with an explicit good rate function depending 
only on the tail distributions of the X i y 

The paper is organised as follows. Section 2 is devoted to the formulation of a large de¬ 
viations principle for the largest eigenvalue x* n of generalised biorthogonal matrix ensembles. 
The examples presented in Section 3 include the random matrix model of disordered bosons 
in [16], biorthogonal Laguerre ensembles considered in [4J as well as Wigner-Dyson ensembles, 
Bogoliubov-de Gennes ensembles and Chiral ensembles. In Section 4 we formulate large devi¬ 
ations principles for a special multiple orthogonal ensemble, the Angelesco ensemble, see [T3 j. 
In Sections 5 we present the proofs of our large deviations principles. 

2. Large deviations for A* of biorthogonal ensembles 

In this section, we will derive a LDP for the bosonic ensemble, where the density of the joint 
distribution of the eigenvalues is of form (11.21) . Obviously, (jl.2[) is a special case of the density 

1 p ( n ) 

g n (Ai,..., A p ( n )) = — |Aj — Aj| |Af — A^| J^w n (Aj) n 1 spw(Ai, • • •, A p ( n )), (2.1) 

n 1 £=1 
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with 0 E N, partition function Z n and continuous weight functions w n : R —)■ Mg". For 0 even, 
£ is a closed subset of [0, oo) while for 0 odd, £ is a closed subset of R. The sequence ( p(n)) n 
must satisfy 

lim ^ = k G (0, oo). 

n—> oo Tl 

Note that for p(n) = n and weight functions ic independent of n, (12.11) subsumes the density 
for the eigenvalue distribution for biorthogonal ensembles as introduced in [3], For 0 = 1, the 
density ( 12 .ip is the same as the density considered in |H] (see formula (4.1) in [9] with 7 = 1 
and (5 = 2). For 0 = 1, p(n ) = n and w(x) = e~ 2 x we also recover the classical GUE. 

Throughout the whole section, we write A f(f) for the set of zeros of a function / : R —» R 
and we assume that the sequence of weight functions (ic n ) n satisfies the following: 

(al) there exists a continuous function w : £ —>■ [ 0 , 00 ) such that 

- #AA(ic) < 00 , Af{w n ) C Af(w) for large n. (al.l) 

— As n —* 00 , w n converges to ic, and log ic n to logic uniformly on compact sets. 
(al. 2 ) 

— logic is Lipschitz on compact sets away from J\f{w). (al.3) 

(a2) If £ is unbounded, then there exists no G N such that 

lim SU p w n (x) = 0 

x —>±00 n>no 

for some fixed e > 0 . 


Moreover we assume that the partition functions satisfy 


where £ is a constant. 


,. 1 , Z n ~ 1 t 

hm - log =: £, 
n-¥ 00 n Z n 


( 2 . 2 ) 


Remark 2.1. We know from [8, Theorem 2.1], that the empirical measure L n = - ^2 5\. of the 

n 3 =1 

eigenvalues of a biorthogonal random matrix satisfying (al) and (a2) obeys a large deviations 
principle with speed n 2 and a good rate function. In some examples the function w n will be of 
the form w n (x) = xntp n (x ) with a fixed a > 0 and some p n (x). But with the first factor xn, 
logic n does not convergence uniformly to some logic on compact sets. This problem, which 
arises in Examples 2.3 and 2.5 in [8], can be fixed easily. We rewrite (12.lj) as 

1 p( n ) 

g n (Ai,..., A p („)) = — ]^[ |Aj — Aj| |A® — A^| A"</9 n (Aj) n l S p(n)(Ai,..., A p ( n )). (2.3) 

n i=l 

n 

Then we obtain that the empirical measure L n = - 22 ■ °f the eigenvalues of this biorthog- 

3 =1 

onal random matrix, where ip n and a limiting function p satisfy (al) and (a2), obeys a large 
deviations principle with speed n 2 and a good rate function. The proof is step by step the proof 
of [8, Theorem 2.1]. 


Remark 2.2. In the proof of our main theorem, one basic step is to integrate out the density 
of all eigenvalues with respect to one and to rewrite some of the remaining parts to a density 
that belongs to an eigenvalue distribution that misses the integrated eigenvalue. Since we do 
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not necessarily have that p(n) = n, we should clarify what we mean by having “one eigenvalue 
less”: We want to analyse q n (X i,..., A p ( n )), which is the joint probability distribution of the 
p(ri) eigenvalues of an n x n matrix. We do so by inserting the probability distribution of p{n) — 1 
eigenvalues. This is not necessarily g n _i, but, since (p(n)) n is a subsequence of the sequence of 
natural numbers (ri) n , there is a p(n — k ) with n > k 6 N and p(n — k) — p(n ) — 1. For practical 
reasons we assume that we have k = 1, i.e. that we have to go back in the sequence (n) n by just 
one step to have one eigenvalue less. This implies that we have p{n — 1) = p(ri) — 1. Further, we 
denote with (with a slight abuse of notation) q n -\{X \,..., A p ( n )_i) the eigenvalue distribution 
of p(n) — 1 eigenvalues. We also have to take care of the normalisation constants of the two 
densities. Their fraction, logarithmised and divided by n, should converge to a constant (see 
(12.2p above). Without restriction, we omit here the first eigenvalue Ai. We also replace, with 
a slight abuse of notation (since, in Equation (12.411 . w p ( n ) and tc p ( n )-i would be more correct), 
w n by w n _ i- This last replacement is technically not necessary; we are at freedom to use any 
transformation of w n that leads to an exponentially equivalent density of the eigenvalues, as 
long as it holds that 


w n (x) 

lim -rr 

"->oo W n -i{X) 


1 . 


(2.4) 


We obtain 


Z n -i 



n i a .-^h a ? 


p(n) 


j=2 


(2.5) 


To be clear: we have some freedom of choice for the functions w n , and, to resume the discussion 
above, it is also possible to use w n -k in the normalisation constant above. Since the sequence 
of the ( w n ) n converges and since we are only interested in the logarithmised and scaled fraction 
of Z n _i and Z n) this does not make any difference. With the same reasoning, the exponent 
n — 1 of w n _i could also be e.g. n — k, with no change in the outcome. 


We will study the asymptotic behaviour of the largest eigenvalue A* : = 
A = (Ai,..., A p ( n )) G The main theorem now reads as follows: 


max^ A j of A for 


Theorem 2.6. Let Ai, ..., \ p ( n ) be the eigenvalues of a biorthogonal ensemble, that is, with 
joint eigenvalue density (12.ip . Under the assumptions formulated above ((al), (a2) and (12.21) ) 
the sequence (A* ) n satisfies a large deviations principle in E with speed n and good rate function 

, = / / (log \x-y\ + log \x 9 - y e \)dp w (y) - logw(x) - C : x > b w 

( oo : x < b w ’ 


where 

C := K J log w(y)dp w (y ) + ^ , 

C, as defined in (12.21) . p w is the limiting measure of the empirical measure of the eigenvalues 
Ai,..., Xp( n ), and b w the right endpoint of its support. Note that ( does not depend on x. 


Note that Theorem 2.1 from [8] ensures that the limiting measure fi w exists. 
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Corollary 2.3. Consider the slightly different density 

1 p( n ) 

g n (Ai,..., A p ( n )) = — ]^[ lAf-Ajfn^)" WAr, ■ ■ ■, A p(n) ), (2.7) 

i=l 

with f3 > 0, 9 e N. Under the assumptions (al), (a2) and (12.2ft i/ie sequence (A*) n satisfies a 
large deviations principle in £ with speed n and good rate function 

_ f -k/3 flog\x 9 - y e \dy w (y)-logw(x) -( ■ x > b w 
' \ oo : x < b w ’ 

where 

C := K J log w(y)dp w (y) + £ , 

£ as defined in H2.2 1) . y, w is the limiting measure of the empirical measure of the eigenvalues 
Ai,..., A p (n), and b w the right endpoint of its support. 


3. Examples 


3.1 Disordered bosons Returning to the bosonic ensemble with density (11.21) . we have as 
weight functions w n (x) = x»e"T =: x « <p n (x) with a fixed «6NU {0}. Remark that r _1 is 
the variance of the independent and normally distributed random variables, that were used to 
construct the stability matrix h for that ensemble, cf. [16] . We have to take the variance r _1 
equal to n^ 1 to be able to obtain a limit of the empirical measures of the eigenvalues. Obviously, 
conditions (al) and (a2) and (12.4)1 are met for <p n {x ) and <p(x) = e~ x , see Remark [2711 We now 
need to verify (12.21) . It is far from trivial to calculate Z n a T for this matrix ensemble. 

For t = 1 and r]i(Xj) = A* -1 , £j(Aj) = Xj 1 and w(Xj) = Xfe~ Xj , we obtain for hj defined 
in the Appendix, that for j > 1 

hjSij = JPi(X)pj(X)w(X) dX = 2 J j\(2j + a)\, 

where the last equality was calculated in [El Equation (5.10)]. Thus applying (16.3p . we get for 
the partition function that 

n— 1 7i—l n 

Zn,a, 1 = C ■ n\ JJ hj = c • n\ 2 J j\(2j + a)! = c ■ 2 n(n_1)/2 JJ j!(2(j - 1) + a)!. 

3=0 3=0 j= 1 

We substitute in the former yj = rXj for all j = 1, ...,n and obtain (with Xj = yj/r and 
dXj = l/rdyf) 



3 =i 


T -3/2n(n-l) T -n(a+l) £ 
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Hence we have proved that 


= C • r -3/2n(n-l)-n(a+l) . 2 n(n-i)/2 + 2 (j _ 1))! . 

j=l 


This gives 


lim — log 1 

n—>oo n Z r , 


= lim - log ( r 3(n-l)+(a+l) . 2-M • 
n—>oo n V 

• 3( n - 1) + (a + 1) 


n!(a + 2 (n — 1))! 


= lim 

n—>oo 


n 


Tl — 1 

log r — lim -log 2 


n—>-oo 77, 


If we choose t = n, we obtain for equation (13.21) 


lim i log 1 = 3(1 — log 2). 


n—>-oo Tl 


(3.1) 


lim — log(n!(o; + 2(n — 1))!). (3.2) 

n—>oo n 


Therefore (j2.2j) is fulhllcd. We may apply Theorem 12.61 on this model. For the rate function, 
we compute the integral 

I(x) = - f (log(z - y) + log(x 2 - y 2 ))deoo(y) - \ogw(x) - (^j log w{y)dg 00 (y) + , 

where g^ is given by (II.3p . With £ = 3 — 3 log2 and ‘'‘ydo^iy) = | the rate function 
becomes 

!{ x ) = ~ J (log(x — y) + log(ic 2 — y 2 )) dg 00 (y) + x — Q + 3-31og2^ . 


3.2 Laguerre biorthogonal ensembles The Laguerre ensembles are a generalisation of the 
biorthogonal ensemble introduced by Luck, Sommers, and Zirnbauer. Take 

, . . . l_ —TX 

S = ( 0 , oo) and w n (x) = x™e ™ , 

with parameter l e No and with r = n and take (12.11) with 0 e N, E = R. and p{n) = n. 
Applying the machinery of biorthogonal polynomials (see (16.3ft in the subsection before and the 
Appendix), we obtain hj = j\9^(9j + l)\ and hence Z n ^i = cn\ ni~o k\9 k [9k + l)\ and 

n— 1 

Z n> Tjt = C ■ r -=^(0+ l)-n(I+l) n | k \Qk(Q k + Z )J . 

k =0 

If the parameter l is constant, we get, with r — n and the use of Stirling’s formula that 

lim — log = Q - 1 - 1 — log 9 . 

n^oo n Z n>n ,i 

This is constant and therefore Assumption 12.21 is fulhllcd. The sequence of weight functions 
w n {\) = \ l / n e~ x =: \ l / n (p(\) converges for n —* oo to w( A) = 93 (A) and, as we have already seen, 
tp fulfils (al) and (a2), see Remark 12.11 Therefore, we can apply Theorem l2.6l to the biorthogonal 
Laguerre ensembles with constant parameter l. That is, we have a large deviations principle 





















KATRIN CREDNER AND PETER EICHELSBACHER 


with speed n for the largest eigenvalue. We omit the calculation of the rate function. If t — n 
and l := l(n ) depends on n such that 

, l(n) t , x 

Inn -=: L E (0, oo) , 

n—»• oo Tb 

where L is a positive constant, the weight function w n converges for n —> oo to X L e~ x =: w( A). 
Here we obtain 

lim — log '^7 L "' / = 6 + L + 1 — log 0 , 

n^oo n ^ n , n ,l 

which is constant and therefore fulfils Assumption 12.21 The assumptions on the weight functions 
are fulfilled as well. Therefore, since all assumptions of Theorem 12.61 are fulfilled, we have a 
large deviations principle with speed n for the largest eigenvalue of the biorthogonal Laguerre 
ensembles. In [4], Borodin mentioned two more prominent classes of random matrices: the 
biorthogonal versions of Jacobi and Hermite ensembles. The calculation for these ensembles 
should be similar to the case we presented, which in turn is again based on the calculations in 

[IE]. 


3.3 The tenfold way Joint densities defined in (12.71) occure is the framework of mesoscopic 
physics, since it subsumes matrix versions of all (ten) classical symmetric spaces, see [9]. It can 
be interpreted as the symmetry classification of disordered fermionic systems. An analogous 
classification for the case of bosons is not completely understood, see, however, the discussion 
in [20, Section 4], For the three classical Wigner-Dyson ensembles we choose f3 > 0, 6 = 1, 
w n (x) = w(x) = e~P x / 4 , p(n) = n and hence n — 1. Obviously assumptions (al) and (a2) are 
fullhlled. Applying Selberg’s integral (see [17, (17.6.7)]), one obtains 



— /3n(n— 1)/4— n/2 n 

n 


i= 1 


r(f) 

r(§) 


With Stirling’s formula we obtain 


lim — log 

n—xx> 77 , 





4' 


Hence assumption (12.21) is fulfilled and we obtain the LDP for the largest eigenvalue, first proved 
in [2J for (3=1. The case (3 = 2 and (3 — 4 are included in [TJ Theorem 2.6.6], but the constant 
£ was not calculated explicitly in [Tj. 


3.3.1. Bogoliubov-de Gennes ensembles. In this section, we prove the large deviations principle 
stated in Theorem 12.61 for four (of five) Bogoliubov-de Gennes (BdG) ensemble (compare with 
Section 4 and the classification table in Section 3 of 0)- We choose 6 = 2, p(n) = n and 
w n (x) = x a / n exp(—-^-A-aX 2 ), where for the four different classes C we choose a = 2, (3 = 2 and 
-0 b = 2 for the class B, a = 0, (3 = 2 and = 2 for the class D, a = 2, (3 = 2 and ifc = 4 

for the class C and a — 1, f3 — 1 and ipci — 4 for the class CL Hence (j2.4[) is fulfilled. We will 
prove that the rate function belonging to Corollary 12.31 for these BdG ensembles is 

/( r) = | P-%1 V t2 - b l dt ■ x > b ™ 

x <b w 


oo 


(3.3) 
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with b w = a j2f) C o 2 fin > 0 being the right endpoint of the support of the limiting measure y w 
of the empirical eigenvalue distribution. This limiting measure is defined as 

^ w{v] = - v 2 ■ ■ (3A 

Note that since 6 is even, E C [0, oo). Therefore, we need to deal only with positive eigenvalues; 
which leads to integrating over only positive values. We obtain I(x) = 0 if and only if x — b w . 
From the upper bound of the LDP and the Borel-Cantclli lemma it follows that P(A* —> b w ^j = 
1 , which is a strong law of large numbers. 


Lemma 3.5. For all BdG ensembles, the scaled ratio of Z n _i and Z n can asymptotically be 
expressed as 


£ 


lim — log 

n—>oo n 




Proof. First, we analyse the partition function Z n . It holds: 


Zn - 






3 = 1 


jpccr' 


r A f)]"II dA 4> 

3 = 1 


where we integrate over the whole space E n C [ 0, oo ] n . Now we use the following transformation 
of Selberg’s integral (cf. [TJ[ Equation (17.6.5)]), which is dehned for positive integers n\ 


OO OO -t 

/ P n n —1 

■■•/ n iA.-A/--nAr'^^=n 

n n l<i<j<n j =1 j =0 


0 0 


F(l + 7 + J7 )F(u + J7) 
T(l + 7) 


(3.6) 


For our purpose of analysing the asymptotics of the partition function, we substitute in (13.61) 

A? 

first A j = ay.j (for some a E M) and then yj = -f, which leads to 


OO OO 


/••■/ n ia?-a?7iiiat”- 1 

^ v, Ki<j<n j=l 

—oo —oo — J — J 




ax _ 7 n(n _i)_„ n ^ n jj £(i + 7 + h)v{v + n) 


3=0 


r(i + 7) 


(3.7) 


Note that the integral in (13.7[) is an even function in all parameters Xj, j = 1,... ,n, so we get 
a combinatorial factor 2" if we integrate over only the positive real axis. We get 


OO OO 


/••■/ n w-AmiAp- 1 

o 0 l<*<3<n 3=1 


- 1 . 


e~i j dXj 


a \-in{n-l)-vn ^ f(l + 7 + j^)T(v + j'y) 


j=0 


r(i + 7) 


(3.8) 
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In order to get a closed formula for the partition function, we take 7 = §, v — § and a = 


This leads to 


Z n ~ 


Z n -l — 


^ a 2 ^n ( „- 1)+S jln . p (1 + . r(g ±l + |(j _ !)) ^ 


n — 1 


i=i 

q + 1 


r(i + f) 


ipcv 2 \ r (i + |j). r («±i + _ 1)) 


n 

t=i 


r(i + f) 


2n 

i>cG 2 


Now we can start the analysis of the behaviour of Z ' 7 1 . We observe 


n—1 


2 \!("-l)(n-2)+S±l(n-l) 


7c ff2 


|n(n-l)+2±ln 


„ N |8(n-l)+ajl ^ f(n-l)(n-2)+2±l(n-l) 


ipcv 2 


n 


(3.9) 


which covers the 11011 -Gamma function expressions of We now consider the part of Z n _ \ /Z n 
that consists of products of Gamma functions. Most of the Gamma functions cancel out each 
other; it remains 

r < 1 + f) (3 10 ) 

r(l + fn)rs±l + f (n - 1 )) 

We apply Stirling’s formula to the Gamma expressions in (13.10p and obtain 

7 4H W y - W 

where the equivalence is logarithmic equivalence , which means for positive numbers (a n ), (b n ) 
that a n ~ h n iff lim^oo ^ log a n / lim^oo - log b n = 1. Note that the numerator of (13.10p is 
logarithmically equivalent to 1. This leads to Equation (13. lOjl being logarithmically equivalent 
to 


(I3T0D ~ 


fin /3(n — l) 

nfj \ 2 f ,Q(n — 1 ) ^ 2 


2 e 

n 


2 e 


2 f3n—/3 

2 f 


2 pn-p Pin- 1 ) 

2\ 2 / n _l\ 2 


) V 2 e 

Now we combine (13.9p and (13.lip and get 


n 


(3.11) 


Z n -1 


n 


'i/jcv 2 


2 2 f P^ccr 


2(3n—f3 

Z\ 2 


2 e 


n 


n 


P(n- 1 ) 


_ X \ f(n-l)b-2)+2±i(n-l)-^ 


This leads directly to 


lim — log “ 1 

n^-oo ft Zrt 


fll I ’P*ca 2 \ , „ , P 
-/Jlogl -2-)+f>+2’ 


where we use the rule of de l’Hospital. 


□ 
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Obviously w(x) 


2 

exp(—^^ 2 ) and hence (al) and (a 2 ) are fulfilled and we obtain 


I(x) 


-/3 / log \x 2 - y 1 1 d/j, w (y) + 


x 


ip c cr 2 


log w(y)dfj, w (y) + £ 


(3.12) 


We define the function 


<h(f, n) := [ log |t 2 - y 2 \dp(y) - - . (3.13) 

J Vc0 2 P 

Consider ^$(t, p w ) for t > b w (hence t 2 — y 2 > 0 for y being in the support of (i w ). For 
convenience, we set c : = 6 2 = 2 i)jc(? 2 f3- 


V~c 

d , . f d , 2 o\ 7 / \ 4 1 f 4 1 2 / -- 4t 

^) = J Tt log (t -y )<Mv) - 7 = J - -■ 

0 


The integral is equal to ^ (2t| — \/t 2 — C7r) and hence ^<F(t, fi w ) = —~ c \/t 2 — c and by the 

X _ 

fundamental theorem of calculus we get 4>(a;, /j w ) = — J -\/t 2 — cdt + $(i/c, y, w ). We calculate 

w c 

4>(i/c, fi w ) = log ( ^ >C 2 — 1- Now it remains to calculate f log w(y) dy, w (y)\ 


j log w(y) d,p w (y) 


4 

C7T 



y 


V’ccr : 


-\Jc — y 2 dy 




2 ' 


Summarising we obtain the desired rate function (|3.3[) . 


3.3.2. The Chiral ensembles. Although the chiral ensembles and the BdG ensembles are closely 
related, there are still a few differences between both models. The most noticeable difference 
is that the empirical measure converges towards a different limiting law. While it converges 
to something resembling a semicircle law in the BdG case, in the chiral case the limiting 
measure of the empirical measure behaves Marcenko-Pastur-like. Another difference lies in the 
number of eigenvalues. While for three out of four of the BdG ensembles, we have n different 
eigenvalues, the chiral ensembles have only s{n ) < n different eigenvalues (specified below). 
Like [9J Section 4], we give the detailed calculations just for the class BDI; the calculations 
for the other two classes AIII and CII are very similar. Consider p(n ) = s(n ) A t{n) and 
9 = 2. Assume, without restriction and for simplicity, that s(n) < t(n ), n G N, and take 

/ > /3(t(n) — s(n))+/3 — 1 ^ 

w n (x) — x~ n~ e 2^r. Denote the partition function by Z n ^ n y As already seen for 

the BdG ensembles, we start with a transformation of the Selberg formula, namely equation 
(13.8ft . We use this equation (with an index shift; we start with j = 1 instead of j = 0) with 
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parameters 7 = /3/2, v — f (t(n) — s(n) + 1 ), a — and n = s(n) and obtain 


J n,s(n) 


°° r s{n) 

— n w-\TniA 

0 0 l<Kj<s(n) t =1 


\P(t(n)-s(n))+p-l e -£z\ 2 jc i\ 


n \ - 2 s ( n )( s ( n )- 1 )-f( t ( n )- a ( Tl )+ 1 ) s ( ri ) 

2^) 

* (n) r (i + If)r(f «n) - «(n) + 1) + (j - i)f ) 


X 


n 


r(i + f) 


2a‘ 


n 


2 s f,(»)(»-»<»)) »M r(i + j|)r(|( n - 2s(n) + j) 


n 

i=i 


r(i + f) 


(3.14) 


To simplify calculations, from now on we set s := s{ri) and t := t(n) (recall that n = t + s). 
We now recall the condition stated in equation (12.41) . that is, lim . = 1. We choose 

n—>■ OO W p(n)-l\ x ) 

s(n — 1 ) = s(n) — 1 = s — 1 and t(n — 1 ) = t(n) — 1 = t — 1 and get 


w n (x) 

^p(n)— 1 (t) 


= \x\ 


, P(t-s)+P-1 P(t-1-S+1)+P-1 

n ti — 2 


-2(,8(t-s)+|8-l) 
— ijq n(n — 2 ) 


This choice of s(n — 1) and t(n — 1) assures that the exponent of x converges to 0, since 
lim s(n)/n — k e (0, 00 ). Consequently, equation (12.4[) is fulfilled. Therefore, we compare s 

n—> oo 

with Z n _ 2 , s _i, where the latter is 


^n—2,s —1 


2a 2 X f ((n- 2 )-( S -l))( S -l) r(1 + J -|)r(§((n - 2) - 2(s - 1) + j)) 


71 — 2 


n 


r(i + f) 


2(T 2 \ |(n-*-i)(>-i) r(i + jf)r(|((n - 2s + j)) 


n — 2 


n 

3 =1 


(3.15) 


r(i + f) 

We now divide Z n - 2 ,s- 1 by Z n)S , using (I3.14j) and (I3.15j) . The first expression in the fraction 

:=B-G (3.16) 


Zn — 2,S — 1 


IS 


B = 


2a 2 \-^n-s)s , 2a 2 \ f (n-s-l)(s-l) f 2(T 2 \^~ n) ( 71 \ 


71 


71 — 2 


71 


2a 2 


Now, we have to deal with the Gamma expressions of and Z n _ 2 , s -i. We get 

r(i + f) 


G = 


r(i + fs)r(|(n-s))' 
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We now apply Stirling’s formula to the Gamma expressions (as for the BdG ensembles) and 
obtain for (13.161) 


Zn—2,s—l 


J n,s 


n \§( n ~ 1 ) f/3s\ aS ff3(n — s)\ 2( " 


2 a 2 


2 e 


n — s /3a 2 


2 \ ~%n 


n 


71 — S 


2e 


-E s 
2 S 


Now, we build the LDP limit and obtain with equation (13.171) 


lim — log Z' 1 2 - b 1 — — ^ ( log(l — k) + k log 


n—>oo n 


J n,s 


K 


1 — K 


n 

2 ^ 


+ log(/3cr 2 ) - 1 


(3.17) 


where we used that lim s/n = k. With lim w n (x) = x^ 1 2l ^e ^ = w(x), the assumptions 

n—>■ oo n—>-oo 

(al) and (a2) are obviously fulfilled. Therefore, we can apply Corollary 12.31 on the largest 
eigenvalue of the chiral ensembles. As rate function, we get 

I(x) =: -<$>(x,n w ) - Ik log w(y)dfM w (y) + f J , 


with ^ 

$(£, n) = /3 k J log |x 2 - y 2 1 dfi(y) + (3{l - 2k) log x ~ ^~ 2 ■ (3.18) 

As in the BdG case, we look at the derivative of $>(x,/j, w ). Note that, for x > b w , it holds that 
x 2 > y 2 because the right endpoint of the support of y w , b w , equals \/b. Therefore, we can omit 
the absolute values in the logarithm: 


d_ 

dx 


$(x,/i w ) = /3 k 


1 


4~ - y 2 ) dfi w (y) - 2L + —— 

dx a z x 


r-Vb 


2x 


<J 2 7T 


V(a 2 -a)(t-!/ 2 )^-4 + /i(1 — ■ 

a z x 


vvw 2 , • (3.19) 

fy/t [x -y )y ° x 

We have y w (y) = - a )(b - y 2 ) with a = 2a 2 ft - yA(l - «)) and 

b = 2cr 2 /3 (§ + — k) j (see 0). We now take a closer look at the integral in Equation 

(13.190 . Standard calculus (or suitable computeralgebra software - we used Mathematica) gives 

2x^(y 2 -a)(b-_ y 2 ))_ = _ tt ^ + ^ _ _ + 

J^e {x 2 - y 2 )y x 


where i = y—1 is the imaginary unit. We used that it holds for the complex logarithm that 
log(— x) = logx + in for all x G M + (this is due to the usual representation of the principal 
value of the complex logarithm and the fact that x has imaginary component 0). Now we 
insert equation (13.20j) in equation (I3.19p and apply the fundamental theorem of calculus on the 
function &(x,y w ) for b w < t < x. Thus, we get 

f x d 

$(z,AG) = / —$(t,ii w )dt + $(b w ,/j, w ) 

J byj 

^ r x i 

=-^ + V( a - t 2 )(b - t 2 ) + t 2 ) dt + $(b w ,y w ), 

Jb w t 
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where <&(b w ,/A w ) is constant. Therefore, we get as a rate function for the BDI ensemble 

I(x) — f -(Vab + yf {a — t 2 )(b — t' 2 ) +t 2 ) dt + c, (3-21) 

° Jb w t 

where c < oo is a constant. Again it follows P(A* —> b w ) = 1. 


4. Large deviations for multiple orthogonal ensembles 


Multiple orthogonal polynomials are a generalisation of orthogonal polynomials in which the 
orthogonality is distributed among a number of orthogonality weights. They appear in random 
matrix theory in the form of special determinantal point processes that are called multiple 
orthogonal polynomial (MOP) ensembles. In (12, H] the appearance of MOP in a variety of 
random matrix models and models related with particles following non-intersecting paths have 
been considered. To a finite number of weight functions n>i,... ,w p on M and a multi-index 
n = (ni,..., n p ) G we associate a monic polynomial Pft of degree n := \n\ := n\ + • • • + n p 
such that 


Pfi(x)x k Wj(x ) dx = 0, for k — 0,..., rij — 1, j — 1,... ,p. 


If Pf, uniquely exists then it is called the multiple orthogonal polynomial (MOP) associated with 
the weights wi,... ,w p and multi-index n. In |T3] the following result was presented. Assume 
that 


det [/) (Tfc)]yfc=i,...,n 

^n 


PI (x k - Xj) 

~l<j<k<n 


(4.1) 


is a probability density function on M n , where the linear span of f \,..., ,/„ is the same as the 
linear span of {x k Wj(x)\k — 0,... , rij — 1, j — 1,... ,p }. Then the MOP exists and is given by 


Pn(x) = E 




4=i 




where the expectation is taken with respect to the p.cl.f (14.111 . which can be interpreted as the 
expectation of the random polynomial n"=i i x ~ x j) with roots x\,... ,x n from a determinantal 
point process on the real line. The p.d.f (14.IK is called a MOP ensemble. It was first observed 
in [3j that random matrix models with an external source lead naturally to MOP ensembles. 
The weights w\,... ,w p are an Angelesco system if there are disjoint intervals Ti, ..., T p C M, 
such that supp(tfj) C Tj, j = 1 ,...,p. In the Angelesco case, det [fj{xk)]j,k=x,...,n is of block 
form and results in 

V / ”-i 

(x k )]j,k=l,...,n = n A^IM) 

i=l A k=1 

with x^ := x Ni _ 1+k G L, AT* = Y?j=i n j (with = 0) and X® = (x[ l \ ... ,Xn}) and 



A(X) = |J (x k - xj) for X = (xi,..., x n ), 

l<j<k<n 
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the Vandermonde determinant. Thus an Angelesco system gives rise to a MOP ensemble, the 
Angelesco ensemble, and the joint p.d.f is 

iWv n A(x">,xw)nn«,(*s»), <4.2) 

i= 1 l<i<7<p i= 1 k= 1 


1 

z„ 


where 

n m 

fc=i j =1 

for X = (x\,..., x n ) and Y = (y 1 ,..., y m ). We now consider the situation that \n\ = n —* oo 
and Uj —> oo for every j = 1 ,..., p in such a way that 

72 ■ 

— —>■ rj for j — 1,..., p (4.3) 

with 0 < rj < 1 and Y2j=i r j — 1- Let us consider varying weights 

Wi(x) = e ~ nVi(a;) (4.4) 


for any i = 1 ,p. Denote by A* := nraxi<fc< nj the the j-th maximal eigenvalue for every 
j — 1,... ,p. We will study the asymptotic behaviour of (Aj,..., A*). 


Theorem 4.1 (LDP for Angelesco ensembles). Assume that every weight function Wi in (j4.4j) 
satisfies assumption (al) and (a2) and assume that (14.3p is fulfilled. Assume moreover assump¬ 
tion ( 12 . 2 p for Z n . Then the sequence (A^,..., A*) n satisfies a LDP on with speed n and good 
rate function 


I(x i,...,x p ) = ^2 r i log \xi-y\ 2 qi*(dy) 

i =1 ^ 

f P 

r i r j / lo s \xi-y\ 1 h*j(y) + TiVi ( Xi ) _ C, 

1 <i<j<P ^ i =1 


(4.5) 


where 


('■= k Y1 / lo § Wi(y)dfj,*(dy) + £. 


2—1 


Here y* = (//*,..., y*) is assumed to be a unique minimiser could in Theorem [ 8 , Section 111] 


Remark 4.2. We would also be able to consider Nikishin ensembles with p > 2 weights, see 
[13] and references therein. This is because the deternrinantal structure of the joint density of 
the eigenvalues [131 (4-14)] consists of Vandermonde-like products. Nikishin interaction arises 
in the asymptotic analysis of eigenvalues of banded Toeplitz matrices as well as in a two-matrix 
model , see [13[ Section 5.4], 


5. Proofs 

This section is devoted to the proof of Theorem 12.61 Theorem 14.11 will not be proved since 
the arguments are very similar, but technically much more involved. 


5.1 Exponential tightness 
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Lemma 5.1. We have 


lim limsup ip n (A* > M) = — oo , 


M—> oo 


n 


For the proof of Lemma 15.11 we need the following technical inequality. 

Lemma 5.2. For all |x| > m&x(b w , 1) and for all A e S, there exists a constant c < oo with 

\x — A| l ^ 61 — A e |w n (A) < c • |x| 0+1 . (5.1) 


Proof. The proof is inspired by the work of mi- First we show that there exist e > 0, T > 0 
and no G N such that for all |A| > T and for all n > n 0 we have 

log<o„(A) < _(l±iKl±£l log(l + A 2 ). (5.2) 

Due to assumption (a2) for unbounded £, there exists an n 0 G N such that 
lim A ^±oo |A|^ +1 ^ K+£ ) sup n > no w n (X) = 0. This implies 

lim ( ( 6> + 1 )( K + ") iQg^ + A 2 ) + log gup Wn (\^ 

A±oo y 2 n>no 

and, since the first summand is positive, it follows that 3 rq G N and 3 T G M such that for all 
n > max(no,ni) and for all |A| > T: 

- log sup (w n (A)) > ( g + 1 )( K + t -) i og (i + A 2 ). (5.3) 

n>n o ^ 

Since the logarithm is a monotonic increasing function, we can cancel the supremum in (15.3[) . 
It follows (15. 2 p for n > max( 77 o,ni) and for |A| > T. Note that this implies that for n > 
inax(no,ni) and for |A| > T that w n (X) < 1. To be complete, we have to cover also the case 
|A| < T. For |A| < T and |x| > 1 we have 

\x — A| \x e — X e \w n (X) < sup 7 n n (A)(|x| + T)(\x\ e + T e ) < c • (xl 61 " 1 " 1 , 

|A|<T 

with c > oo being a constant. Now we consider the left-hand side of inequality (15.ip and start 
with a case distinction: For |x| > |A|, we have 

|x — A| \x e — A 0 |u; n (A) = \x — X\\x e — A 0 | exp (logiu n (A)) 

< \x — A| |a ; 61 — A e | exp ^ + + ^ log(l + A 2 )^ (5.4) 

< \x — Alla;' 9 — A 0 | < 4|x||a: 0 | < 4|x| 0+1 , 

where we used (I5.2p for the first inequality in (15.4p and for the second inequality that ( e+ 1 h K+£ ) > 
0, log(l + A 2 ) > 0. The third inequality is due to the fact that \x e — A 0 | < \x e \ + |A 0 | < 2\x e \ 
for |x| > | A|. The last inequality uses 0 G N. For |x| < |A|, we start again with inequality (I5.4j) 
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and obtain: 


\x 


\x e — X e \w n (X) < \x — A| \ x® — A^| exp ( ——————-log(l + A 2 ( 


< 4|A|‘ l+1 exp(- (9+1) i' I + ' ;) log(l + V 


— 4exp ( (6 + 1) log |A — — — + ^ !og(l + A 2 ] 


(5.5) 


Now, since |A| > |x| > max(fe^, 1) and since c := (6,+1 A h+ b > 0, the term —clog(l + A 2 ) 
dominates the exponent in (j5.5[) . Therefore, the whole exponent is bounded by 1 and we obtain 
fl53p < 4 - 1 < A\x\ e+1 . □ 


Proof. (Proof of Lemma 15. 1 1) 

We take a closer look at the following density: 


Pn(A * n >M) = 


n P( n ) P(. n ) 

/ / IWV" n |Ai-A,||A»-A"in<iA 

M 3 =1 1 <i<3<P{n) j = 1 


M £p(n)-l 

3 = 1 



OO 

r r 


p(n) 

pin) 

/ / 

W n (Xl) 

” II l Al - AIlN - A ?l • 

n 

M J>(n)-1 


3 = 1 

3 = 2 

p{n) 




JJ w n 

-i(A :) n 

- 1 n ^-•mia? 

- A 1 

3= 2 





(.O n ( A y 


p(n) 

n 

i =2 


Z n -1 


M sp (^)- 1 


°n P( n ) 

/ / ^AirniA.-vw-A 

J J j =2 


p(n) 


* If rfP »-.( A 2 ■ ■ ■ . 

j = 2 w n-l\*j) 


(5.6) 


where we added ^=2 Z" l(V)^- 1 ^ ie s ^ e P anc ^ substituted P„_i in the second step. Now, 


write the last product in 

(|5.6j) as 


n »„(a t Y ft ( 

w n {Xj) N 

n- 1 P(n) 

\ TT 

j=2 W n-l(^j) n ~ 1 

Wn—l(^j ) y 

I 11 

3=2 
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Therefore, we get for equation (j5.6[) that 


p„(a; > m) 



(5.7) 


< c 


Z n -i 


ir n (A 1 ) n |Ai| (0+1)(p(n)_1) 


M SPt™)- 1 
p(n) 

><n 

3 =2 


w n(Aj) 

l^n— 1 (Aj 


n—1 


c ®’n-i(^2 • • • Ap( n ))dAi, (5.8) 


where we used Lemma [5.21 for the inequality with constant c = L Now, we take a closer 
look at a part of the integrand of (15.8j) : 


^n(Ar) n |A 


n \ \ |(0+l)(p(n)-l) _ 


^(AOIA,!^)^^ 


|Ai 


i—(0 +i)(e+^:) 


For large n, the fraction is near its limiting value n. For big |Ai|, we know from assumption 
(a2) that the first bracket of equation (15.91) is small. Therefore for large n and large |Ai| we 


have: (j5.9[) < |Ai| n ( 0+1 )( e +i). Since we assumed in equation (I2.4[) that lim u Wn ^^ = 1, we 


have lim J {'' _ 


p(n) ( WnjXj) 


j ^ V — 1 {^j ) 


n —1 


= 1. Therefore, (15. 8 j) becomes 


flEED ^ c- 


= c • 


Z n -1 
Z n 

Z n ~\ 


| Ai |-n(0 +1 )( e +i) dXi / dP n _ l( A 2 . . . A p (n)) 


M 

oo 


Vp(n)-i 


| Ai |-n(«+l)(s+i) (jAi 


(5.9) 


M 


for large values of n and |Ai|. The last step is due to the fact that P n -i is a probability. Now 
we integrate (15.9p out with respect to Ai and get 

Z„_ i 1 


flElD = c- 


-M l ~ n {d+ l ){e+l) _ 


Z n 1 — n(9 + l)(e + -) 

Now we take the LDP-limit and obtain 

1 1 i Z i 

lim sup — log P n ( A* > M) < lim sup — log c + lim sup — log —^— 

n —»oo U n—>oo Tl n —>oo Tl Z/ n 

1 1 

+ hm sup - log -- 1 

ih oo u 1 — n{6 + 1)(£ + —) 

1 


lim sup — log M 1_n(e+1)(£+ ») 

n—¥ oo Tl 
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The first summand is zero, the second equals £ due to assumption (12.2[) . The third summand 
has size O(^p) and thus converges to 0 for n —* oo, and the last summand is 

1 . //iiiw , i\ 1 — 

limsup — logM 1 n ( 6,+1 )(£+J = limsup-— logM = —clogM 

n—> oo ^ n—> oo ^ 


for a positive constant c. Therefore, since £ is a constant, we obtain the result. 


□ 


Since the density of the eigenvalues behaves symmetrically for x —» oo and x —> —oo, we can 
analogously prove the following lemma: 

Lemma 5.10. Define A^ := min^i A j. Then we have 

lim limsup — logP ri (A^ < —M) = — oo . 

M—tco n _^ QO Tl 


5.2 / is a good rate function 
Lemma 5.11. I(x) is a good rate function. 


Proof. Since we assumed in (al) that w(x) is continuous, I(x) is continuous in x on the interval 
(b w , oo) and lower semicontinuous in x = b w . For any x, y G R it holds log \x — y\ < log(|x| + 
1) + log(| 2 /| + 1). Hence for any probability measure // 

J log \x — y\ + log \x e - y e \dfji{y) 

< f log(|x| + 1) + \og(\y\ + 1) + log(|x e | + 1) + log(|/| + 1) dy(y) 

= f logOI + l)(|x 0 | + 1) dfi(y) + j log(| 2 /| + 1)(|/| + 1) dfi(y) 

= log(|x| + l)(|x 6> | + 1) + J log ( 12/ 1 + 1)(|/| + 1) dy(y) . (5.12) 


With (15.12j) it follows for the rate function 

I(x) = — k J log \x — y\ + log \x 6 — y e \ dy{y) — logw(x) — c 

> -«log(|x| + l)(|x"| + 1) - f./log(|i,| + l)(|y»| + l)dM»)-log»W-c. 

(5.13) 

Again, we neglect those terms that are constant concerning x (i.e. the integral and c in (15 . 13H ). 
We also omit the absolute values, since we are interested in large positive values of x. We apply 
inequality (15.2)1 . which gives information about the behaviour of w(x) for x — * oo, and obtain 
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for x > rnax{fc,„, 1} that 


— k log(x + 1) (x e + 1) 


logiu(x) > 

> 


— K log(x + 1) (x 6 + 1) + + + l 0 g(l + x 2) 

— K log X^ 0+1 ^ + log ^ + + ^ log(l + X 2 ) 


= O (log(a;- K(m) )) + £> (x (0+1) ( K+£) ) 

= 0 (logx e(e+1) ) . 

This converges to oo for x —» oo since e > 0 and 9 G N. Therefore, J(x) is strictly increasing 
on [ b w , oo), it has compact level sets and is thus a good rate function. □ 


5.3 Proof of a weak LDP In this section, we prove a weak LDP for A* (for a definition see 
[6J Section 1.2]). Since / is a good rate function, it suffices to show that for any x < b w , 

limsup — logP n (A* < x) = —oo (5.14) 

n—>-oo Tl 

and for any x > b w , 

lim — logP n (A* > x) — — I(x) , (5.15) 

n—>• oo Tl 

where b w is the right endpoint of the support of the limiting measure. The reason we need 
to prove (I5.14p and (15.151) lies in the following consideration: [6j Theorem 4.1.11] states that 
for the proof of a weak LDP, we do not need to show the LDP upper/lower bound for every 
compact /open set belonging to the space E. Instead, it suffices to look at a base of the topology 
ofE C R. We choose as a basis the set of open and closed intervals A := {[x, y\, ( x , y) : x < y}. 
Still, we need to verify 


I(x) = sup 

{ AeA-.xeA} 


lim sup — log P n (A* G Al) , 

n—>• oo Tl 


(5.16) 


which is, for closed sets A = [x, y] and for all x < y , equivalent to 

— inf I(t) = limsup — logP n (A* G [x,y]) . 

te[x,y] n—>■ oo Tl 

But this is easy to see with equations (I5.14j) . (15. 15ft and the fact that / is increasing on [ b w , oo): 
We briefly distinguish three cases: if x < y < b w , then limsup A logP n (A* G [x,y]) = —oo due 

n—>• oo 

to equation (15. 14ft . If x < b w < y, we use the fact that in this interval lies the “typical” value 
of A*, therefore limsup 4-logP n (A* G [x, y]) = 0. Finally, if b w < x < y, we get with equation 

n—>oo 

(I5.15j) that lim sup 4 logP n (A* G [x, y]) — — inf^^ ^ J(t). With the continuity of /, we get that 

n^f- oo 

we have the same limits if we apply the same calculations to open intervals (x, y ) instead of 
closed intervals [x, y\. Therefore, we obtain a weak LDP for A* as soon as we have proved 
equations (j5. 14[) and (I5.15p . which is what we do in the rest of this section. 

n 

We first prove (15.14p . We know from [51 Theorem 2.1], that the empirical measure L n = - ^2 5\. 

n j =i 

of the eigenvalues of a biorthogonal random matrix obeys a large deviations principle with speed 
n 2 and a good rate function. Since I is a good rate function, it achieves its minimum. Now we 
fix x < b w . Since A* < x < b w (this means that there is no eigenvalue of X n in the interval 
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(x, b w ]), we find a bounded continuous function / G C&(£) with f f dL n = 0 but f f dp, w > 0. 
Therefore, we conclude with P, Theorem 2.1] that the probability of {A* < x } is exponentially 
decaying (with speed n 2 ) and therefore it follows equation (I5.14p . We split up the proof of 
Equation (j5.15[) into two steps. We first show that —I(x) is an upper bound of the scaled 
probability on the left-hand-side of (j5.15f) . then we prove that —I{x) is also a lower bound. 


5.3.1. Proof of the upper bound. Observe that we have, for any M > x > b w , the following 
inequality: 


p„(a; > x) < p n (A; g [x, m}) + p n (A; > m) . ( 5 . 17 ) 

Since we proved that large values of A* are exponentially negligible, we choose M large enough, 
i.e. that the first term in the right hand side of inequality (15. 17(1 is exponentially small, and 
also that the same holds for the minimal eigenvalue A^, see Lemma [5.101 Hence we only need 
to deal with the probability that A* G [x, M ]. Further it holds 


k / log \x-y\+ log \x° - y \d/i w (y) + logw(a;) + k log w(y)dn w (y) 


> sup (a / log |t y\ T log I r - y°\dii w (y) + log w(t) + k \ log w(y)dn w (y) ) . 

fE[M, oo) 


(5.18) 


which we will apply later on in the proof. Since the eigenvalues \j, 1 < j < p(n), are exchange¬ 
able, it holds P n (A* G [x,M]) < nP n (Ai G [x,M\, |Aj| < MVj > 2). Now we consider 


Pn(Ai G [x,M], |Aj-| < MMj > 2) 



pin) 

A j | 'ui^(Aj) d\ 2 * * * dd\\ 

j =i 



xdP n -i(A2, • • •, A p(n)) d\\ , 


(5.19) 


with Im '■= (S H (— M, M]) p ^ 1 , and where we replaced P„_i as in equation (15.6p . Now we 
define, for t G and fi supported on [— M,M ], the function 


n 


p(n ) — 1 




log \t-y\ + log 1 1 9 - y e \dfi(y) + 


log w n (t) + 


log w n (y)dn(y). (5.20) 
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We consider the first part of the integrand in equation (15. 19[) : 


p(n) 

™4Ai) n nK(Ai)|Ai-A j ||A;-A?|] 

3 =2 

p(n) p(n) 

= exp [n log w n (Xi) + '^logw n (X j ) + ^[log|Ai - Xj\ + log | Af - Aj|] } 

3 =2 3 =2 


= exp < (p(n) — 1) 


log |Ai -y | + log |A L -y \ dL n _ x {y) + 

+ J log w n (Xj) dL n _ x {y) 

= exp{(p(n) - l)$„(Ai,L n _i)}, 


n 


p(n) — 1 


log w n { Ai) 


(5.21) 


p(n) 

where i := —bj- Y 5^. is the empirical measure of the p(n) — 1 eigenvalues A 2 ,..., X p ( n )- In 
3 =2 


p(n) 

the third line, we used the identities Y l°g | Af — A®| = (p(n) — 1) J log | Af — y 9 \ dL n _ x {y) and 

3 =2 


P{n) 

Y logu; n (A j) = {p{n) — 1) f log w n (y) dL n _ x {y). Therefore, with equations (I5.20p and f!5.2ip . 
3 =2 

we obtain for Equation (15.19P 


P„(A 1 G[x,M],|A J |<MVj>2) 


Z n -1 


M 



3 (p(n)-l)<I' n (Ai,L rl _i) 


Z n -1 


x I M 
M 


p(n) 

n 

3 =2 


) 

Wpi— 1 (^j) 


n—1 


d!Pn-l(A2, • • • ? ^p(n)) 6^1 



(p(fi)-l)^ n (Ai,L n _i)+(n-l) £ lQ gf 

; i= a V dF n _ 1 (X 2 ,...,\ p(n) )dX 1 


X I M 

M 


Zn- 1 f f (p(n)-l)^ ra (Ai,L n _i)+(n-l) 



s Im 


dP n -l(A2, • ■ ■ j Ap(n)) dAi 


M 

= / J e Wd- 1 )^(h,d (i (p ri _ 1 oL-i 1 )( /i ) ( iA 1 , (5.22) 

1 Mi(En[-iW,Af]) 


where 


$n(f,A*) 


<h n (f, p) + (n - !) J log p(dy). 


(5.23) 
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Now, in equation 05.220 . we split the domain of integration of the inner integral. 

Mi(En[-M,M]) = [Mi(En[-M,M])nB(fj, w ,8)] 

U[M!(S n [-M, M]) n (B(fi w ,5)) c ] 
= B m (ji w , 5) U [Af!(S n [-M, M}) n (£(//„, 5)) c ]. 

This leads to 05.221) being equal to 
Pn(A]_ G [x,M],\Xj\ < MWj > 2) 


Z n -i 


M 

/ J 

M 




e (p(n)-l)*„(Ai>/0 rf(p n _ 1 o 




^n-1 


(/i + / 2 ). 


(5.24) 


We first discuss the integral J 2 , which (as we shall see) converges to — oo on an LDP scale. To 
show this, we prove an upper bound for <f> and integrate over a larger set than that appearing 
in the definition of J 2 . We choose Ai G [x, M] and /j G .Mi(E fl [— M,M ]). Then we have 
log |Ai — Xj\ < log(2M) and log | Af — A^| < log(2 M 6 ). Therefore, for as defined in 05.201) . 

$n(Ai,)U) < log(2M) + log(2M e ) + —--- sup logw n (t) + sup logw n (t), 

P\ n ) -l- te[x,M] te[-M,M\ 

and therefore we have for <f> n , as defined in 05.231) . 


<MAi,/i) < log(2M) + log(2M e ) + 


n 


sup log w n (t) 
p[n) — 1 te[x,M] 


+ sup log w n (t) + (n - 1 ) sup log lL ' 1 ^ . 
te[-M,M] te[-M,M] w n -i (t) 

Thus, we have an upper bound for <f> n that is independent of A 2 ,..., A p ( n ) and that we can place 
before the inner integral. Further, we have for the domain of integration of the inner integral 
of J 2 that {jMi(S f) [— M, M]) fl (B(/j, w ,S)) c } C {(B, 5)) c }. This leads to 


1 2 < exp <| log 2M + log 2M e H———-- sup log w n {t) + sup logw; n (£) 

p(n) — 1 t&[x,M] te[-M,M] 


n 


u i Wn (*) 

+ (n - 1) sup log- — 


—i ° L n — i)(h ^ B([i w , 5)). 


Now it follows that 

/} n 

J 2 = 2M ■ 2M° sup ta n (t) p(n )- 1 sup w n (t) sup 


tE[x,M] 


Wn(t) 


(n-1) 


^ \ / / \ 
te[-M,M] te[-M,M] \ w n-l\t). 

x(P n _! O L~\)(n (£ B(fi w ,6)) 


< 4M ( ' 6,+1 ^ 1+ ^ K+£ ' ) pG"-i ) +2 • e “ c (<5)n 2 


( 5 . 25 ) 
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for some M,n large enough and for some positive constant c(<5). The second step is due to the 
fact that we exchanged the logarithm and the supremum for w n . This is allowed for n large 
enough because we are on a compactly supported interval where w n converges uniformly to 
the continuous function w. The third step uses assumption (a2), the boundedness of [x, M] 
and [— M,M] as well as the fact that (P„_i o obeys a LDP with speed n 2 , see |8] 

Theorem 2.1]. Since the order of the first factor in (j5.25[) is smaller than e n2 , the estimation is 
complete because now J 2 vanishes on the LDP scale. Remember that the LDP for A*, which 
we want to prove, has speed n. So we only have to deal with the first integral of equation 
(I5.24p and, of course, the factor Z n _ x jZ n , applying [BJ Lemma 1.2.15]. As an aside, we also 
get with 0 Lemma 1.2.15] that it suffices to consider P n (A* £ [x,M]) for the calculation of 
lirnsup A log P n (A* > x), cf. inequality (15.17jl and the fact that in the LDP scale, P n (A* > M) 

n—> oo 

converges to — 00 . Therefore we have: 


lim sup — log P n (A* > x) 

n—>-oo Tl 


= lim sup - log P n (A* £ [x, M ]) 

n—> oo Tl 

< lim sup — log (1 1 ) 

n—>-oo Tl Zj n 


< £ + lim sup — log / 1 . 

n—> OO Tl 


(5.26) 


The last line is due to assumption (12.21) . Now we shall take a closer look at the integral 
I\. It holds for a probability measure p, £ Adi(S) with p £ B M (fi w ,S) that <E> n (Ai, L„_i) < 
supfieB M 0^,9) ^n(Ai ,p), which is independent of A 2 , ■ ■ ■, A p ( n ). Hence we obtain 


M 

h < y’ e (p W- 1)sup A^M(M.,«)^DnA) dAl J d(P n _! o 

% 5^) 

M 

< J e (p(n)_1) SU P*"(Ai.A) 


< exp{(p(?z) - 1) sup $„(f, n)} ■ (M - x ); (5.27) 

te[x,M\ 


the second step follows because P n _i °L~\ is a probability measure. In the third step we build 
the supremum over all allowed values of Ai, extract the integrand (which is now independent 
of the integration variable) out of the integral and integrate out over 1. Now we look at I\ on 
the LDP scale. It follows with inequality (15.270 : 


lim sup — log /1 

n —>-00 Tl 


< 


lim sup 

n—> 00 


(pH -1) 

n 


sup $ n (t, /i) + lim sup — log(M 

B M (f! w ,S) n-¥ 00 n 

tP[x,M] 


x) 


= hi sup lirnsup $ n (Ai, fx). (5.28) 

neB M (tM»,S) n-¥ 00 
\i£[x,M] 


Now we step back and consider for n —> 00 . On compact intervals, logu7 n converges 
uniformly to log re due to assumption (al). Therefore, if we insert this in the definition of 
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fl5.23f) . we get 

limsup <i» n (Ai, n) = \og\Xx~ y\ + log y e \d/j,(y) + log w(X 1 ) 
n—>-00 J 

Wn(y) 


log w(y)dy(y) + limsup / log 


w n -\(y) 


n— 1 


dy(y) , 


where the last summand converges to 0 for n —> 00. We define as limit of 4> n and <f> n the new 
function 


<h(Ai,/i)=: j log |Ai - y\ + log I A® - y\dn(y) + k Mogiu(Ai) + J log w(y)dfi(y) (5.29) 
and obtain, if we insert this in (15.281) . 

lim sup — log I\ < k sup <h(Ai, y). 

n-+o o n / ieB M (n w ,6) 

Ai £[x,M] 


By (15.261) . we have 


lim sup — log P n (A* > x) < £ + lim sup — log I\. 

n—>■ 00 ^ n^oo Tl 


and therefore 


lim sup — log P n (A* > x) < C, + k lim sup <E»(Ai, //). 

n -¥00 U y,£B M {iMn, 5 ) 

Ai e[x,M] 


(5.30) 


Note that the calculation for the scaled and logarithmised integral I 2 is also independent of 
6 , since the calculations there hold for all 6 > 0. Since the measures /i that are inserted in 
come from a probability space endowed with a topology that is compatible with weak 
convergence, we have 


lim sup <f>(Ai,/i) = sup $(Ai,/!„,)= sup <h(Ai,/i w ). 

AiS[a;,M] AiSp.oo) 

AiS[x,M] 


(5.31) 


The last equality is due to (15.181) . We now look at the relation between the function <f> and 
the rate function / from Theorem 12.61 ft holds for all t > b w $(t,[jb w ) = — A _1 (/(t) +0- Since 
we proved that I(x) is strictly increasing on [b w , 00) and continuous on (b w , 00), we have that: 
<E>(x, y w ) is continuous for x > b w (and lower semicontinuous for x = b w ) and <3>(x, y w ) is strictly 
decreasing for x > b w . Therefore, the supremum in equation (15.311) is attained in the leftmost 
value that Ai can attain, namely, x. That implies, inserted in 05. 301) . for all x > b w 

lim sup - log P n (A* > x) < i + k$(x, y w ). 

n^-oo Tl 

This is the desired upper bound we wanted to prove. 


5.4 Proof of the lower bound In the last step of the proof of Theorem 12.61 we prove that our 
rate function —I is a lower bound for the scaled probability on the left-hand-side of equation 
(IS.lSp . We fix y > x > r > b w and 5 > 0. Since {Ai G [x, y], niax^ |Aj| < r} C {A* G [x, M]} 
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for large M and y close to x, we look at the following probability: 


P n (Ai e [x,y\, max |Aj| < r) = 


3 =2 



x J jMi(Sn[— r,r\) 




(compare with (15.2211 ). For Ai G [x,y], the integrand has the following lower bound: 


(5.32) 


exp ((p(n) - l)<h n (Ai, n)) > exp 


{pin)- 1) inf <1 n (t,fj) 
te[x,y) 


This leads to a lower bound for (15.32)1 : 

(0211 > [ V d\ 1 f e W n )- 1 '> i ° f *ei*.v]*''( t ’rid(W n _ 1 o L-\)(n) 

J x J _Mi(£n[—r,r]) 

^n JB r (fJ,w,S) 

> —^—(y — x) exp | (p(n) — 1) inf $ n (t, /r) ip n _ 1 (L n _i e 5 r (^, 5)), 

Ai l t£[x,y\ J 

flEBr {flw ,S) 


where in the second step we built the second integral over the smaller set {L n _i E B r (y w , 5)} C 
{L n _i E Afi(Sn [—r, r])} and integrated out the first integral. In the third step, we once more 
made the integrand independent of the integration variables and put it outside the integral. 
Now, we build the LDP limit: 


lirninf — logP n (Ai E [x, y], max | A,-1 < r) 

n—> oo Tl j =2 


> lim inf — log ” 1 + lim inf — log(r/ 


x 


n—>■ oo 77 , 


n—> oo Tl 


+ lim inf 

n—>■ oo 


p(n) — 1 


n 


inf $„(*,//) 

t&[x,y\ 
flE.Br (flw ,< 5 ) 


+ liminf-logP n _i(L n _i E B r (fj, w ,6)) 

n—>-oo Tl 

= £ + lim inf ——-- inf $ n (t, /i), 

n-> oo n 

flE.Br (flw 5 ^) 


(5.33) 

(5.34) 


where we used assumption (12.2[1 and the fact that the second and the last summand of inequality 
(I5.33jl are zero. Observe that for the last term this follows from the large deviations result 
for the empirical measure that provides P n _i(L n _i T- B(y w ,5)) — » 0, therefore P n _i(L n _i E 
B(/j, w ,S)) —* 1 and it holds logP n _i(L n _i E B(y w ,5)) 0. We just look at the second 
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summand of (j5.34[) and fill in the definition of <f> n as defined in equation (15.231) : 
lim inf ~ 1 inf $ n (t,/j) 

n-> oo n te[x,y\ 

(f-lw j< 5 ) 

T)('fl) — 1 // 77 / 

= liminf —- inf / log \t - y\ + log \t e - y e \ dfi(y) + —- -- log w n (t) 

n—too n te[x,y\ \J p[n) — 1 

+ I log w n (y)dn(y) + (n - 1) J log d ^(y)^j 

>n inf ( f log\t~y\ +log \t e - y e \dfi(y) 
te[x, y } \J 

y€B r (Hyj ,< 5 ) 

+a _ 1 logiu(i) + Jlo g w(y)d M ) (5.35) 

— k inf $(£, /i). 

te[x,y\ 


For inequality (I5.35F we used again (as in the proof of the upper bound) the uniform conver¬ 
gence of log w n to login on compact sets as assumed in assumption (al) and exchanged the 
Limes inferior and the inhmum. The last equality is just the definition of <f>, see equation f|5.29p . 
Therefore we have 

1 p(ri') 1 

lim inf — logPAAi G [x, y ], max |A 7 j < r) > £ + k; lim inf — log( inf $(t, /i)). 

n—)■ oo U j =2 n^-oo Jl te[x,y\ 

fl^Br j^) 

Now we again build the limit 8 \ 0 on both sides, which does not affect all terms but one, and 
observe that the function <f> is continuous on [x,y] x A / fi([—r, r]) for y > x > r > b w . We also 
let y \ x and obtain 

lim inf — logP n (Ai G [x, y ], max I A,j < r) > k<L(x, /i w ) + £ , 

n—>oo n 3 =2 

which is the desired lower bound. We therefore proved Theorem 12.61 


6. Appendix 


Borodin [4j starts with a set of n real random variables Ai,..., A n with values in the interval 
/CM and with joint probability density function 


Qn ( Al, • 


• •; A n ) 


det(?7i(A :; ))j ! j = i i ... !n det(^(Aj))j J= i,... in 

X/TJ, 


• n-M ’ 

i=i 


( 6 . 1 ) 


where is the normalisation constant and y and £ are real-valued functions. We suppose 
that we are able to biorthogonalise the families {r]i} and {£*} regarding the weight function w, 
i.e. that we have an inner product (£,?/) = f T £(A)y(A)w(A) dX. This entails that we find two 
other families of functions, {Pi(A)}j l =1 and {q , j(A)}” =1 where the p^ and g* are monic (i.e. the 
leading term has prefactor +1) polynomials in A of degree i — 1, and that are orthogonal (not 














28 


KATRIN CREDNER AND PETER EICHELSBACHER 


necessarily orthonormal) with respect to the weight function, i.e. we have pi £ Span{£i,..., £ n }, 
qi £ Span{? 7 i,..., rj n } and 


hjSij := ( Pi,qj) = I Pi(X)q j (X)w(X) d\ Vi,j . 

For the construction of such polynomials, see e.g. m Section V.A], where the authors used a 
modified Gram-Schmidt algorithm to obtain the desired polynomials. Define now the matrix g 
with 


9i,j '■= j rji{X)£j(X)w(X) d\ , 

and assume it to be not singular. Now we have everything at hand to formulate two propositions 
that help to calculate the normalising constant of biorthogonal ensembles, see PI Proposition 
5.5] and 0 Equation (24)]. 


Proposition 6.2. For the normalisation constant of general biorthogonal ensembles with den¬ 
sity (16.Ill , we have 

71—1 

Z n = nl det g = c ■ n\ h 3 , (6.3) 

3 =0 


where c is a constant. 
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